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Quantum information science attempts to ex- 
ploit capabilities from the quantum realm to ac- 
complish tasks that are otherwise impossible in 
the classical domain [lj. Although sufficient con- 
ditions have been formulated for the physical re- 
sources required to achieve quantum computa- 
tion and communication [2], there is an evolv- 
ing understanding of the power of quantum mea- 
surement combined with conditional evolution for 
accomplishing diverse tasks in quantum informa- 
tion science [3, 4, 5j. In this regard, a significant 
advance is the invention of a protocol by Duan, 
Lukin, Cirac, and Zoller (DLCZ) 0| for the re- 
alization of scalable long distance quantum com- 
munication and the distribution of entanglement 
over quantum networks. Here, we report the first 
enabling step in the realization of the protocol of 
DLCZ, namely the observation of quantum cor- 
relations for photon pairs generated in the col- 
lective emission from an atomic ensemble. The 
nonclassical character of the fields is evidenced 
by the violation of a Cauchy-Schwarz inequality, 
9i2($t) ^ <?i.i52,25 with Qij as normalized corre- 
lation functions for the two fields = (1,2). 
As compared to prior investigations of nonclas- 
sical correlations for photon pairs produced in 
atomic cascades [7] and in parametric down con- 
version 0, our experiment is distinct in that the 
correlated (1,2) photons are separated by a pro- 
grammable time interval St, with St ~ 400nsec in 
our initial experiments. 

The theoretical proposal of DLCZ 6] is a probabilis- 
tic scheme based upon the entanglement of atomic en- 
sembles via detection events of single photons in which 
the sources are intrinsically indistinguishable, and gener- 
ates entanglement over long distances via a quantum re- 
peater architecture 0. The DLCZ scheme, with built-in 
quantum memory and entanglement purification, is well 
within the reach of current experiments and accomplishes 
the same objectives as previous more complex protocols 
that require still unattainable capabilities |9lllCj|. 

In our experiment, we demonstrate a basic primitive 
integral to the DLCZ scheme. Specifically, an initial 
write pulse of (classical) light is employed to create a 
state of collective atomic excitation as heralded by pho- 
toelectric detection of a first photon 1. After a pro- 
grammable delay St, a subsequent read pulse interrogates 
the atomic sample, leading to the emission of a second 
(delayed) photon 2. The manifestly quantum (or non- 



classical) character of the correlations between the initial 
write photon 1 and the subsequent read photon 2 is veri- 
fied by way of the observed violation of a Cauchy-Schwarz 
inequality for coincidence detection of the (1,2) fields 0. 
Explicitly, we find \gl a {5t) = (5.45 ± 0.11)] £ [91,152,2 = 
(2.97 ±0.08)], where are normalized correlation func- 
tions for the fields and St = 405 ns is the time sepa- 
ration between the (1, 2) emissions. The capabilities real- 
ized in our experiment provide an important initial step 
toward the implementation of the full protocol of DLCZ, 
which would enable the distribution and storage of en- 
tanglement among atomic ensembles distributed over a 
quantum network. Extensions of these capabilities could 
facilitate scalable long-distance quantum communication 
6] and quantum state engineering [Tl| . For example, 
by employing spin-polarized samples in optical-dipole or 
magnetic traps [l^. it should be possible to extend the 
interval St to times of several seconds. 

Our experiment arises within the context of prior work 
on spin squeezing 0, Q] , and in particular on atomic 
ensembles where significant progress has been made in 
the development of methods to exploit collective en- 
hancement of atom-photon interactions provided by op- 
tically thick atomic samples [Tl 111 Irl 111 111 |2C|. 
Instead of homodyne or heterodyne detection of light 
as used in spin-squeezing experiments |l8l HI l20l ]. 
the DLCZ scheme involves photon counting techniques, 
which present stringent requirements for broad band- 
width detection and for the suppression of stray light 
from the atomic ensemble. 

As illustrated in Fig. ^ an optically thick sample 
of three-level atoms in a lambda-configuration is ex- 
ploited to produce correlated photons via the following 
sequence. With atoms initially prepared in state \a) by 
optical pumping, a laser pulse from the write beam tuned 
near the \a) — ► \e) transition illuminates the sample and 
induces spontaneous Raman scattering to the initially 
empty level \b) via the \e) — » \b) transition at time t^ x \ 
The write pulse is made sufficiently weak so that the 
probability to scatter one Raman photon into the pre- 
ferred forward propagating mode ip^(r, t) is much less 
than unity for each pulse. Detection of a photon in the 
mode ipW(r, t) produced by the |e) — > \b) transition re- 
sults in a single excitation in the symmetrized atomic 
level I b). In the ideal case, this coherently symmetrized 
state is p 
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Although the initial detection of photon 1 generated by 
the write beam is probabilistic, the detection of photon 1 
results in the conditional state \$>ia) with one collective 
atomic "excitation" . This excitation can subsequently be 
converted into an excitation of the light field with high 
probability "on demand" with a specified emission direc- 
tion and a programmable pulse shape 0, 0, 0, |2^| ■ In 
order to achieve the conversion from atoms to field, a 
laser pulse from the read beam tuned near the \b) — > |e) 
transition illuminates the atomic sample, thereby affect- 
ing the transfer \b) — > \a) for the sample with the accom- 
panying emission of a second Raman photon 2 on the 
\e) —> \a) transition. For an optically thick atomic sam- 
ple, photon 2 is emitted with high probability into a spec- 
ified mode V (2) (r, t) offset in time by i (2) = + St. The 
spatial and temporal structure of the modes tp^ 1,2 \r,t) 
are set by the geometry of the atomic sample and by the 
shape and timing of the write and read beams |21| . In 
our experiment, the modes of the (write, read) beams are 
spatially mode matched, with measured visibility greater 
than 95% for the case of equal frequency and polariza- 
tion. The time delay St is limited in principle only by 
the coherence time between the levels |o) and \b), which 
is long in practice. 

The atomic sample for our experiment is provided 
by Cesium atoms in a magneto-optical trap fMOT) [l2l| . 
where the Cs hyperfine manifolds {\6Si/2,F = 
i),\6Sy 2 ,F = 3),|6P 3/2 ,F' = 4)} correspond to the 
levels {\a), \b), |e)}, respectively. As illustrated by the 
timing diagram in Fig. ^ the MOT is chopped from 
ON to OFF with At = 4/js. In each cycle there is a 
"dark" period of duration 1/js when all light responsible 
for trapping and cooling is gated OFF, with less than 
0.1% of atoms measured to remain in the F = 3 level 
at this stage. The j th trial of the protocol for single 
photon generation is initiated by a write pulse which is 
resonant with the 6Si/ 2 ,F = 4 — > 6Ps/2,F' — 4 tran- 
sition at frequency 1^4,4 and that has duration ~ 51 ns 
(FWHM). A critical parameter for the experiment is the 
resonant optical thickness 74.4 of the atomic sample [2l]| . 
We measure 74,4 ~ 4 — 5 for cw excitation, corresponding 
to an attenuation of intensity exp(— 74,4) in propagation 
through the MOT. 

The write pulse generates forward-scattered (anti- 
Stokes) Raman light around frequency W3.4 from the 
F' = 4 excited level to the F = 3 ground level (|e) — > \b)) 
that is directed onto a single-photon detector Dl. Af- 
ter a variable delay St, the read pulse illuminates the 
sample, with this pulse tuned to the d>Si/2,F — 3 — > 
6P 3 /2iF' = 4 transition at frequency 0^4 with duration 
~ 34ns (FWHM). Raman (Stokes) light generated by the 
read pulse around frequency W44 from F' = 4 to F = 4 
(|e) — > I a)) is directed onto a second single-photon detec- 
tor D2. 

By interchanging the frequencies for optical pumping 
of the filter cells described in Fig. ^ the (write, read) 
beams can be detected at (Dl, D2) in place of the (1, 2) 
fields. An example of the resulting pulse profiles accumu- 



lated over many trials {j} is presented in Fig. [21 where 
the origin in time is set to coincide with the approximate 
center of the write pulse, with the read pulse following af- 
ter a delay ~ 415ns determined by external control logic. 

With the filter cells set to transmit the (1, 2) photons 
to the (Dl, D2) detectors, respectively, we record his- 
tograms of the numbers (ni(t),ri2(t)) of photoelectric 
events versus time, which are also displayed in Fig. [21 
For the data presented here, the intensity of the write 
pulse is kept low (~ 10 3 photons per pulse), resulting in 
a time lag for the onset of the n\ (t) counts in Fig. [21 As 
discussed in the Supplemental Information, the probabil- 
ity polite t° generate an anti-Stokes photon 1 within the 
solid angle of our imaging system is p^l ite — 10~ 2 per 
pulse. 

The read pulse is about 100 times more intense than 

the write pulse, leading to high efficiency £3-»4 ~ 0.6 

(2) 

for the transfer of population \b) — > \a), with p read — 

(3^iPwlite f° r the Stokes photon 2. Examples of the re- 
sulting detection events ^(i) are shown in Fig. [21 In 
contrast to the behavior of n\(t), the intense read beam 
generates 712(2) counts promptly. More extensive investi- 
gations of the timing characteristics of the emitted fields 
(1,2) will be part of our subsequent investigations, in- 
cluding the relationship to electromagnetically induced 
transparency (EIT) |1[2|. 

A virtue of the DLCZ protocol is its insensitivity to 
a variety of loss mechanisms, including inefficiencies in 
transport and detection of the (1,2) photons. However, 
in an actual experiment, various non-ideal characteristics 
of the atom-field interaction (as in our MOT) do lead to 
deterioration of correlation for the (1,2) photons (e.g., 
imperfect filtering and/or background fluorescence as de- 
scribed in the caption of Fig. [21 and in the Supplemen- 
tal Information). Fortunately there exists a well-defined 
border between the classical and quantum domains for 
the (1,2) fields that can be operationally accessed via 
coincidence detection, as was first demonstrated in the 
pioneering work of Clauser ■ 

As illustrated in Fig. \%)p, electronic pulses from detec- 
tors (Dl, D2) are separately gated with windows of du- 
ration T = 60ns centered on times (t^ , t^ ) correspond- 
ing to the approximate peaks of the (n\(t) , n2(t)) pulses 
shown in Fig. [21 Photoelectric events that fall within 
the gate windows are directed to a time-interval analyzer 
(TIA) configured in a standard fashion for measurement 
of photoelectric correlations j2(| . For a start event from 

Dl within the interval vp ±T/2 for the j th trial of the ex- 
periment, the TIA records the times of stop events from 

(2) 

D2 within successive intervals t£' ± T/2. Over many 
repetitions of the experiment, we thereby acquire time- 
resolved coincidences ni^r) between the (1,2) fields, 
both within the same trial k = j and for subsequent 
trials k = j + l,j + 2... (i.e., a start event from trial 
j around time tj and a stop event from trial k around 
time tk, where tk = tj + (k — j)At for k = j, j + 1, . . .). 
By a 50-50% beam splitter, the field 1 can be directed to 
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detectors (Dl, D2), and then in turn the field 2 to (Dl, 
D2). We thus also acquire the time-resolved coincidences 
ni,i(r) and n 2 ^ 2 (r). 

Figure [3] displays an example of data accumulated in 
this manner for coincidences n a ^(r) between the (1,2), 
(1,1), and (2,2) beams, with successive peaks sepa- 
rated by the time between trials At = 4/is. Note that 
there is an excess of coincidence counts in each of the 
initial peaks for joint detections from the same trial 
(t < At) as compared to n a ^(r) from different tri- 
als (r > At). This excess is shown more clearly in 
the plots in the right column, which expand the time 
axis from the left column in Fig. |3| Here, data from 
successive trials k = j + 1, ...j + 10 have been off- 
set to t < At and then averaged for comparison with 
n aji a(r) from the same trial j by introducing the quan- 
tity m a ,p(r) = J2kt)°+i n a ,p{r + (k - j)At). As dis- 
cussed in the Supplemental Information, statistical inde- 
pendence for trials with k =/= j is enforced by the exper- 
imental protocol of reapplying the MOT and repuming 
beams after each trial. 

From the data in Fig. |3| we determine the total num- 
ber of coincidences N a .p = J2{ Ti } n a,p( T i) with (a, (3) — 
(1. 2) obtained by summing over time bins {t/\ for detec- 
tion within the same trial j, and M a ^ = X){r fc } m a,p( T k) 
obtained from start and stop events from different trials 
(j ^ k). Fields for which the Glauber-Sudarshan phase- 
space function is well-behaved (i.e., classical fields) are 
constrained by a Cauchy-Schwarz inequality for the var- 
ious coincidence counts [Supplemental Information and 
Ref. namely 

[h,2(5t)f < 51,152,2 , (2) 
, ~ _ Ni,i ~ _ W22 ~ / rjA _ JVi, 2 

where g hl = , 52,2 = ^ , g h2 (6t) = 

For the data displayed in Fig. [3J we find 51,1 = 
(1.739 ± 0.020) and 32,2 = (1-710 ± 0.015), in correspon- 
dence to the expectation that the (1, 2) fields should each 
exhibit Gaussian statistics with 5^1 = 92,2 = 2 for the 
protocol of DLCZ in the ideal case, but here degraded by 
diverse sources of background counts (see Supplemental 
Information). By contrast, for the cross-correlations of 
the (1,2) fields, we record gi, 2 (St) = (2.335 ± 0.014), 
with 5t = 405ns. Hence the inequality of Eq. for 
classical fields is strongly violated, namely [g^ 2 (St) = 
5.45 ± 0.11] ^ [51, 152,2 = 2.97 ± 0.08], where all errors 
indicate the statistical uncertainties. This violation of 
the Cauchy-Schwarz inequality clearly demonstrates the 
nonclassical character of the correlations between pho- 
tons (1,2) generated by the (write, read) beams. More- 
over, as discussed in more detail in the Supplemental 
Information, the measured coincidence rates in Fig. [31 



explicitly document the cooperative nature of the emis- 
sion process. Overall, we estimate that the probability 
pi 9 ' for coincidence of the (1,2) photons due to collec- 
tive atomic excitation as described by the state \§ia) is 

roughly pi 9 ^ ~ 10~ 4 for each trial j, referenced to the 
output of the MOT. 

The temporal extent of the photon wave packet ip(r, t) 
for the (1,2) photons is also of some interest. To in- 
vestigate this issue, we have carried out the experiment 
with expanded gate windows of duration T — 140ns that 
then encompass the entire domains over which counts 
rii(t) and n 2 (t) are observed in Fig. [21 In this case, 
we record g 1A = (1.72 ± 0.04), 52,2 = (1-52 ± 0.05), and 
g 1>2 (6t) = (2.45 ± 0.10), now with 5t set to be 320ns. 
The classical inequality of Eq. is once again not satis- 
fied; [gl 2 (5t) = 6.00 ± 0.50] ^ [51,132,2 = 2.61 ± 0.11] . 
These results with T = 140ns also confirm that dead- 
time effects do not play a significant role in the current 
experiment. 

As described in the Supplemental Information, the vi- 
olation of the Cauchy-Schwarz inequality of Eq. [21 in 
the ideal case can be much larger than we have ob- 
served, namely [gi, 2 (St)) 2 /[51, 152,2] - [(1 + P)/(^P)f > 
1, where p <C 1 is the excitation probability. In our ex- 
periment, the size of the violation of the inequality was 
limited mostly by uncorrelated fluorescence from indi- 
vidual atoms in the atomic sample. This contribution 
will be made smaller in future experiments by moving to 
off-resonant excitation, which necessitates higher optical 
density. There is also a significant limitation due the 
presence of the leakage light from the read pulse. This 
classical pulse is only 9GHz away from the single-photon 
field 2 of interest, and is filtered by a factor exceeding 
10~ 9 . To achieve even stronger violation of the inequal- 
ity, we must further improve the filtering capability. 

Our observations of nonclassical correlations between 
the (1, 2) photons represent the first important step in the 
realization of the protocol of DLCZ [|| for scalable quan- 
tum communication with atomic ensembles, although it 
is not yet sufficient for realization of the full protocol. 
Beyond the nonclassical correlations, our experiment also 
demonstrates successful filtering of the various fields and 
collective enhancement by the atomic ensemble, all of 
which are critical for realization of the full quantum re- 
peater protocol. More generally, the capabilities that 
we have demonstrated should help to enable other ad- 
vances in the field of quantum information, including 
implementation of quantum memory |22l l27j ] and fully 
controllable single-photon sources [28| . which, combined 
together, help to pave the avenue for realization of uni- 
versal quantum computation Q . 
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FIG. 1: A simplified schematic of the experiment is presented, 
with (a) providing a diagram of the apparatus and (b) giving 
the timing sequence for data acquisition. Further details are 
as follows, (a) Write and read pulses propagate sequentially 
into a cloud of cold Cs atoms (MOT) , generating pairs of cor- 
related output photons (1,2), with controlled separation St. 
Fields with frequency near that of the \a) <-> |e) ([&) «-» |e)) 
transition are colored red (blue) here and in the subsequent 
two figures. The write and read pulses have orthogonal po- 
larizations, are combined into a single input at PBS1 (PBS - 
polarizing beam splitter), and are then focussed into the Cs 
MOT with a waist of approximately 30/im. The output fields 
are split by PBS2, which also serves as the first stage of filter- 
ing the (write, read) beams from the (1, 2) fields. For example, 
field 2 is transmitted by PBS2 to be subsequently registered 
by detector D2, while the read pulse itself is reflected by 90° 
at PBS2 and then blocked by an acousto-optical modulator 
that serves as Gate 1. Further filtering is achieved by passing 
each of the outputs from PBS2 through separate frequency 
filters each of which consists of a glass cell of Cs vapor op- 
tically pum ped to place atoms into either 6Si/2,F = 3 or 
F = 4 23]. The small residual reflected (transmitted) light 
of the write (read) pulse from PBS2 at frequency u)4, 4(0^3,4) 
passes through a filter cell with atoms in the F — 4(3) level. 
It is thereby strongly attenuated (> 10 6 ), while the accom- 
panying Raman-scattered light as photons 1(2) at frequency 
0)3,4 (0)4,4) is transmitted with high efficiency (~ 80%). Trans- 
mission efficiencies from the MOT to detectors (Dl, D2) are 
both about 30% for light with the spatial shape of the write 
and read beams and of the correct polarization. (Dl, D2) 
have overall quantum efficiencies of approximately 50% (pho- 
ton in to TTL pulse out), (b) Gating windows for the joint 
detection of photons (1,2) are centered at times (tj\tj) for 
the j th trial of the experiment during intervals when the MOT 
is OFF. 
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FIG. 2: Normalized singles counts rii(t) are shown for the 
write, read, and (1,2) fields. The pulses around t = are 
from detector Dl for the write beam n w (t) (solid trace) and 
for photon 1, n\(t) (points). The pulses around t = 410ns 
are from detector D2 for the read beam n r (t) (solid trace) 
and for photon 2, nzit) (points). Note that in addition to 
the symmetrized excitation, each write pulse also transfers 
several hundred atoms into the F = 3 level due to spontaneous 
emission from its near-resonant character. However, atoms 
transferred into F — 3 via spontaneous decay are spatially 
uncorrelated, so that their contribution to the signal from the 
read channel is strongly suppressed (by roughly the fractional 
solid angle collected, Sfl/iir ~ 4 x 10~ 5 ) as compared to the 
signal from single-atom excitations of the form 0- 
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FIG. 3: Time-resolved coincidences n a ,^(r) between the 
(1, 1), (2, 2), and (1, 2) fields are displayed versus time delay r. 
Left column - n Q ,/3(r) is shown over 11 successive repetitions 
of the experiment. Right column - The time axis is expanded 
to a total duration of 250ns with r = set to the center of the 
gating window (t^,^,^) for (nn, ri22, niz), respectively. 
The larger peak n a> p(r) corresponds to detection pairs from 
the same trial j, while the smaller peak m a ,s(r) is for pairs 
from different trails as defined in the text. Typical acquisition 
parameters are as follows. Detectors (Dl, D2) have average 
count rates of about (400/s, 250/s), respectively, while back- 
ground counts with no MOT present are about 100/s. Counts 
due to the MOT itself (with write and read beams blocked) 
are less than (10/s, 20/s) for (Dl, D2). Dark counts with 
the inputs to the fibers blocked are less than 5/s. All these 
numbers are for the gated-output mode of data acquisition as 
in Fig. El with T = 60ns. 



